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LONG-TERM  GOALS 

The  long  term  goal  of  this  project  is  to  determine  optimal  sampling  strategies  for  drifting  buoys,  in 
order  to  enhance  prediction  of  particle  motion  in  the  ocean,  with  potential  applications  to  ecological, 
search  and  rescue,  floating  mine  problems,  and  design  of  observing  systems. 

OBJECTIVES 

To  develop  methods  for  the  objectives  of  ODDAS  DRI,  which  is  to  preserve  the  sampling  volume  of 
an  acoustic  sensor  array,  and  replenish  when  gaps  in  the  observational  volume  occur  due  to  the 
dispersion  of  the  sensors  with  the  underlying  Eulerian  flow  field. 

APPROACH 

The  work  is  based  primarily  on  stochastic  and  dynamical  models  of  particle  motion  and  data 
assimilation  strategies.  It  also  involves  the  use  of  coastal  circulation  models  and  processing  of 
oceanic  data  such  as  drifter  positions,  ocean  surface  currents,  and  wind  field. 

WORK  COMPLETED 

1)  Publication  of  the  Lagrangian  Analysis  and  Predictability  of  Coastal  and  Oceanic  Dynamics 
(LAPCOD)  book  by  the  Cambridge  University  Press. 

2)  Publication  of  three  papers  (Caglar  et  ah,  2006;  Haza  et  ah,  2007a;  Chin  et  ah,  2007) 
for  which  most  of  the  research  was  carried  out  during  the  previous  grant  year,  but  their 
revision  still  required  some  effort  and  time. 

3)  Analysis  of  targeted  drifter  trajectories  from  the  winter  cruise  of  the  DART  (Dynamics  of  the 
Adriatic  in  Real  Time)  observational  program,  and  publication 

of  a  paper  (Haza  et  ah,  2007b). 

4)  Analysis  of  relative  dispersion  from  synthetic  drifter  trajectories  derived  from 
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the  Navy  Coastal  Ocean  Model  (NCOM)  configured  in  the  Adriatic  Sea,  and 
submission  of  a  manuscript. 

RESULTS: 

1)  Results  from  targeted  drifter  launch  strategies  during  DART : 

Pis  (Griffa,  Ozgokmen)  and  postdoc  (Angelique  Haza)  have  participated  in  the  drifter  release 
component  of  Dynamics  of  the  Adriatic  in  Real  Time  (DART)  experiment  conducted  jointly  by 
NATO-NURC  (led  by  Michel  Rixen)  and  NRL  (led  by  Jeff  Book).  DART  involved  two  trials,  namely 
DART06A  in  March  2006  and  DART06B  in  August  2006.  We  focus  on  observations  during 
DART06A  here,  concentrating  on  a  set  of  surface  drifters  launched  during  a  high  resolution 
hydrographic  ship  survey.  Drifter  launches  were  specifically  directed  based  on  maxima  of  finite-size 
Lyapunov  exponents  calculated  from  the  NCOM  configured  in  the  Adriatic  Sea  to  provide  two-day 
hind  and  forecasts.  The  specific  goal  was  to  choose  positions  that  would  maximize  the  particle 
spreading  and  therefore  the  spatial  coverage  of  sampling.  To  the  knowledge  of  the  authors,  this  is  the 
first  time  that  drifter  launches  during  an  observational  program  have  been  guided  in  near  real  time 
based  on  modem  coastal  model  output  and  Lagrangian  techniques.  A  main  objective  of  the  multi- 
institutional  DART  experiment  in  the  coastal  area  of  the  western  central  Adriatic  is  to  study 
mesoscale  instabilities  arising  in  the  Western  Adriatic  Current  (WAC)  near  the  Gargano  Cape.  A 
suite  of  different  measurements  were  used  in  conjunction  with  a  real-time  modeling  effort.  As  part  of 
DART06A,  a  total  of  12  surface  drifters  (CODE  and  SVP  types)  were  successfully  launched  between 
March  1 1  and  23,  2006. 

The  main  sampling  goal  for  the  drifters  was  to  achieve  a  good  coverage  of  the  DART  area.  To  this 
end,  launches  were  planned  in  the  two  dynamically  distinct  areas  namely  in  the  WAC  and  along  the 
inter-gyre  boundary  offshore  from  the  Gargano  Cape.  The  exact  launch  locations  and  times  were 
constrained  by  logistical  considerations  including  weather  conditions  and  the  pre-planned  ship  track. 
The  launch  locations  in  the  WAC  were  determined  a-priori  based  on  historical  information  and 
required  no  real-time,  model-based  direction  effort.  Three  directed  launches  were  made  targeting  the 
region  off-shore  from  the  Cape.  The  resulting  drifter  trajectories  provide  a  good  coverage  of  the 
region,  satisfying  the  general  goal  of  the  sampling  strategy.  This  is  in  part  due  to  directing  drifter 
launches  along  the  out-flowing  branch  of  identifiable  Lagrangian  boundaries.  This  strategy  optimizes 
the  relative  dispersion  of  drifters,  the  overall  data  coverage,  and  the  sampling  of  high  kinetic  energy 
features  in  the  flow  field.  A  prerequisite  for  implementing  this  strategy  is  a  methodical  identification 
of  the  Lagrangian  features. 

A  practical  method  to  identify  high-dispersion  regions  and  mixing  boundaries  upon  which  to  base 
drifter  launch  strategies  is  the  computation  of  the  Lagrangian  spatial  structures  produced  by  the  local, 
finite-size  Lyapunov  exponents  (FSLE).  Model  FSLE  fields  were  calculated  from  sets  of  5  synthetic 
trajectories  centered  at  every  grid  point  in  the  NCOM  model.  The  FSLE  map  of  the  DART  region  at  a 
particular  day  is  obtained  by  advecting  a  total  of  1 14,250  particles  both  forward  in  time  from  the  2 
day-forecast  velocities,  and  backward  in  time  from  the  2  day-hind-cast  velocities.  Once  the  location  of 
a  strongly  hyperbolic  region  was  determined  from  isolated  intersections  of  the  FSLE  fields  from  two- 
day  NCOM  forecasts  and  hind-casts,  a  subset  of  10-20  synthetic  drifters  were  launched  in  this 
selected  area  and  advected  by  the  2-day-forecast  fields  as  a  preliminary  test.  Finally,  2-3  drifter 
launch  locations  were  chosen  to  straddle  the  in-flowing  structure  and  to  lie  as  close  to  the  out-flowing 
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pathway  as  the  ship-track  would  allow.  These  locations  were  then  communicated  to  the  ship  on  a 
real-time  basis. 


Figure  1:  The  initial  2-day  trajectories  for  real  (gray:  upper  panels,  green  and  purple: 
lower  panels)  drifters  launched  on  March  16,  March  19  and  March  23  (launch  position 
indicated  by  circles).  Superimposed  are  the  FSLE  computed  at  launch  time  (a)-(c),  and 
synthetic  drifters  released  in  regular  arrays  (d)-(f).  Red  (blue)  dots  indicate  initial 
(final)  positions  of  the  synthetic  drifters. 

Launch  set  1:  March  16,  buoys  14715  and  44924 : 

Due  in  part  to  the  low-intensity  winds,  a  well-defined  hyperbolic  point  (HP)  is  visible  from  the  FSLE 
map  (Fig.  la).  Two  locations  are  selected  along  the  track  corresponding  to  the  out-flowing  Lagrangian 
branches  from  the  HP,  based  on  the  FSLE  extrema.  As  a  result,  the  drifters  diverged  immediately  and 
propagated  in  opposing  directions.  The  success  of  this  launch  is  due  in  part  to  the  nearly-stationary 
position  of  the  HP  over  the  following  two  days,  and  effectively  real  time  implementation  of  the 
suggested  launch  locations.  The  synthetic  trajectories  also  confirm  the  polarized  nature  of  the 
Lagrangian  dispersion  in  the  vicinity  of  the  launch,  as  all  synthetic  particles  are  transported  along  only 
two  out- flowing  branches  (Fig.  Id). 

Launch  set  2:  March  19,  buoys  44925  and  4492: 

Launch  locations  were  initially  chosen  based  on  the  model-forecasted  FSLE  fields  for  March  17,  to 
straddle  an  in-flowing  branch  of  an  identified  Lagrangian  structure  boundary  near  a  HP  with  the 
expectations  of  rapid  divergence  of  the  drifter  pair.  But,  logistical  problems  delayed  the  actual  drifter 
launch  until  March  19.  As  shown  in  Fig.  lb,  the  model  accurately  predicts  the  rapid  reaction  of  the 
flow  to  the  changing  wind  conditions  and  the  movement  of  strong  intersections  of  FSLE  maxima 


3 


away  from  the  previously  determined  launch  locations.  As  such,  the  actual  drifter  launch  occurred  in  a 
region  predicted  by  the  model  to  have  minimal  relative  dispersion. 

Launch  set  3:  March  23,  buoys  44928  and  44930: 

In  this  case,  the  deployment  time  was  constrained  by  the  finite  extent  of  the  observational  program, 
with  the  latest  possible  launch  date  being  March  24-25.  While  the  model  did  not  predict  the  presence 
of  a  strong,  isolated,  hyperbolic  point  along  available  ship  tracks  on  March  23,  Fig.  lc  indicates  the 
presence  of  a  Lagrangian  structure  boundary  in  the  launch  region  as  given  by  distinct  maxima  in  the 
backward-in-time  FSLE  (the  red  curve).  The  actual  drifter  deployment  placed  one  drifter  in  the  interior 
of  this  structure  and  the  other  outside.  This  deployment  leads  to  advection  in  nearly  opposite 
directions  and  a  rapid  separation  of  the  observed  trajectories  in  keeping  with  the  overall  predictions  of 
the  model  (Fig.  If).  Drifter  44930  followed  the  pronounced  northward  pathway  outside  the  model 
FSLE  structure  with  remarkable  accuracy  (Fig.  3c).  Initially  the  other  drifter,  44928,  propagated 
southwards  with  the  structure  as  predicted  by  the  model.  The  observed  drifter  track  shows  appreciable 
slowing  before  further  southern  propagation  due  to  the  break  up  of  the  transport  barrier  blocking 
southward  pathways  one  day  after  the  launch  (not  seen  in  launch  time  FSLE  plot).  This  change  in 
structure  and  opening  of  advective  pathways  is  presumably  due  to  the  influence  of  the  Bora  wind 
event. 

We  emphasize  the  remarkable  spatial  complexity  and  rapid  change  of  the  model  FSLE  fields  (Fig.  la- 
c),  forced  and  maintained  mainly  by  winds  and  mesoscale  turbulent  coherent  flow  structures.  While 
the  envelope  of  synthetic  drifter  trajectories  and  the  observed  trajectories  indicate  significant 
differences,  the  model-derived  FSLE  fields  (resulting  from  averages  over  both  time  and  ensembles  of 
trajectories)  are  apparently  robust  enough  for  accurate  short-time  prediction  of  the  observations. 

To  our  knowledge,  this  is  the  first  time  that  targeted  drifter  launch  strategies  have  been  employed  in 
real-time  during  an  observational  program.  Overall,  we  conclude  that  the  accuracy  of  the  surface 
Eulerian  fields  from  state-of-the-art  operational  coastal  models  appears  to  have  reached  a  level  that  is 
adequate  enough  for  model-based  Lagrangian  diagnostics  to  be  used  in  real-time-directed  drifter 
launches  in  observational  programs. 

2)  Estimates  of  relative  dispersion  in  the  Adriatic  Sea: 

Synthetic  drifter  trajectories  computed  from  velocity  data  produced  by  a  high-resolution  NCOM 
model  are  used  to  investigate  the  scaling  of  relative  dispersion  and  the  distribution  of  Finite-Scale 
Lyapunov  Exponent  fields  in  the  Adriatic  Sea.  Relative  dispersion  quantifies  how  fast  particle  pairs 
separate  from  one  another,  and  thus  it  is  a  key  metric  for  the  purposes  of  ODD  AS  DRI.  The  effects  of 
varying  degrees  of  spatial  and  temporal  filtering  of  the  input  Eulerian  velocity  fields  on  the 
Lagrangian  statistics  are  investigated  in  order  to  assess  the  sensitivity  of  such  statistics  to  model  error. 
It  is  shown  that  the  relative  dispersion  in  the  model  Adriatic  circulation  is  generally  super-diffusive, 
scaling  nearly  ballistically  in  close  agreement  with  Lagrangian  observations  from  a  limited  set  of 
drifters. 
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Figure  2:  Spatial  distribution  of  FSLEs  forward  six  days  in  time  starting  from 
September  6,  2002  (positive)  and  backward  six  days  in  time  (negative)  calculated  from 
(a)  the  raw  field,  (b)  with  1.5-day  time  filtering,  (c)  with  5-day  time  filtering,  (d)  with 
30-days  time  filtering,  (e)  with  spatial  smoothing  filter  scale 
of  2.5  km,  (f)  5  km,  and  (g)  10  km. 
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It  is  found  that  the  large-scale  dispersion  is  dominated  by  persistent  separation  regions  and  the 
controlling  influence  of  the  Western  Adriatic  Current.  Temporal  filtering  with  averaging  windows  up 
to  monthly  time  scales  only  affects  the  relative  dispersion  at  scales  smaller  than  20  km  without  altering 
the  overall  scaling  regime.  In  contrast,  spatial  smoothing  at  scales  as  small  as  5  km  significantly 
reduces  relative  dispersion  at  all  scales  up  to  100  km.  While  no  clear  exponential  regime  is  observed  in 
the  full  data  set,  a  distinct  flattening  of  the  FSLE  curves  are  observed  for  larger-scale  spatial 
smoothing.  The  role  of  chaotic  advection  in  determining  the  small-scale  relative  dispersion  in  such 
cases  is  examined. 


Figure  3:  Effect  of  space-smoothing  on  (a)  the  relative  dispersion,  and  (b)  on  the  FSLE 
curve  for  six  different  values  of  spatial  filtering  parameter.  The  trajectories  near  the 
coastline  have  been  discarded  to  avoid  spurious  dispersion,  and  the  difference  in 
relative  dispersion  and  FSLE  is  highlighted  by  the  new  (black)  and  original  (gray)  full 

field  curves. 
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IMPACT/APPLICATIONS 


The  investigation  of  the  predictability  of  particle  motion  is  an  important  area  of  study,  with  a  number 
of  potential  practical  applications  at  very  different  scales,  including  searching  for  persons  or  valuable 
objects  lost  at  sea,  tracking  floating  mines,  ecological  problems  such  as  the  spreading  of  pollutants  or 
fish  larvae,  design  of  observing  systems  and  navigation  algorithms. 
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